New Class of Eigenstates in Generic Hamiltonian Systems 
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In mixed systems, besides regular and chaotic states, there are states supported by the chaotic 
region mainly living in the vicinity of the hierarchy of regular islands. We show that the fraction 
of these hierarchical states scales as h~ a and relate the exponent a = 1 — I/7 to the decay of 
the classical staying probability P(t) ~ t~"' . This is numerically confirmed for the kicked rotor by 
studying the influence of hierarchical states on eigenfunction and level statistics. 



Typical Hamiltonian systems are neither integrable nor 
ergodic jl| but have a mixed phase space, where regular 
and chaotic regions coexist. The regular regions are or- 
ganized in a hierarchical way j|] (see, e.g., Fig. 2a) and 
chaotic dynamics is clearly distinct from the dynamics of 
fully chaotic systems. In particular, chaotic trajectories 
are trapped in the vicinity of the hierarchy of regular is- 
lands. The most prominent quantity reflecting this, is 
the probability P(t) to be trapped longer than a time t, 
which decays as 
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P(t)~r\ 7 >1 
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in contrast to the typically exponential decay in fully 
chaotic systems. While the power-law decay is universal, 
the exponent 7 is system and parameter dependent. The 
origin of the algebraic decay are partial transport barri- 
ers j|, e.g., Cantori, leading to a hierarchical structure of 
the chaotic region 0. Quantum mechanically the clas- 
sical algebraic decay of P(t) is mimicked at most until 
Heisenberg time ||. 

Even after two decades of studying quantum chaos, the 
search for quantum signatures of this universal power-law 
trapping is still in its infancy: In fact, only conductance 
fluctuations of open systems have been investigated so 
far. It was semiclassically derived that these fluctua- 
tions should have a fractal dimension D — 2 — 7/2 0, 
which was confirmed in gold nanowires [h, semiconduc- 
tor nanostructures |]] , and numerics [ |To[ | . Quite recently, 
a second type of conductance fluctuations in mixed sys- 
tems has been discovered numerically 11,1^], namely iso- 
lated resonances. There the classical exponent 7 seems 
to appear in the scaling of the variance of conductance 
increments, surprisingly, on scales below the mean level 
spacing, what is not understood so far. Thus, even for 
this subject there is a lack of basic understanding. 




FIG. 1. Husimi representation of a regular, hierarchical, and 
chaotic eigenstate of the kicked rotor (K = 2.5, h = 27r/1000) with 
solid lines showing KAM-tori of the classical phase space. The 
lower part shows the distribution p(r]) of the average density 7] an 
eigenstate has in the chaotic region far away from the islands (green 
shaded area corresponding to 25% of sites in g-representation) . The 
mean of the Gaussian fit determines l// c h, which is clearly distinct 
from l/p ch . 

In this paper we present consequences of the classi- 
cal P(t) ~ t f in closed quantum systems, namely the 
appearance of a new class of eigenstates. Different from 
the well studied regular and chaotic states, they are sup- 
ported by the chaotic region, but predominantly live in 
the vicinity of the regular islands |p^| . As for decreas- 
ing h they move deeper into the hierarchical structure of 
the chaotic region, we like to call them hierarchical states. 
We show that the fraction /hi er of such states scales as 
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which is confirmed numerically for the kicked rotor by 
studying their influence on eigenfunction and level statis- 
tics. We deduce relation (0) by combining the finite reso- 
lution of quantum mechanics for a given fi with the sim- 
plest model J|| describing classical transport in a mixed 
phase space. We are encouraged by the success of this 
approach, and think that it will help in understanding 
the quantum signatures of a mixed phase space. 

We motivate and numerically verify our results using 
the well known kicked rotor, which is a paradigm for a 
generic Hamiltonian system |Jl4| and has a time evolution 
described by the map: 

q n +i = q n + p n 

Pn+i = Pn + K sing n+ i . (3) 

Its quantum properties are determined by the time evo- 
lution operator for one period fl5| , |l6| 

U = exp f-^Kcosq\ exp • ( 4 ) 

Fig. 2a shows a typical classical phase space, which is 
governed by one big chaotic invariant set with fractional 
phase space volume p cn and islands of regular motion 
of total fraction p xeg jL7| arranged in a hierarchical way. 
Quantum mechanically, however, we find three types of 
eigenstates of U (Fig. 1) : There are 'regular' states living 
on KAM-tori of the regular islands and there are 'chaotic' 
states extending across most of the chaotic region as first 
described by Percival p8| . In addition, there are eigen- 
states supported by the chaotic region, that are of a third 
type. They predominantly live in the vicinity of the reg- 
ular islands with only a small contribution in the main 
part of the chaotic sea. These states are separated from 
the main chaotic sea by partial barriers of the classical 
phase space, e.g., Cantori or stable and unstable mani- 
folds ||] . For decreasing h these states live deeper in the 
hierarchy of classical phase space and thus we like to call 
them hierarchical states. Denoting the fraction of the 
three types of states by / rcg , / c h and /hier, respectively, 
one has 

/rcg + /ch + /hier = 1 = Prog + Pch ■ (5) 

Obviously for any given Ti not all regular islands are re- 
solved and therefore f Te g(K) < Prog holds. We determined 
the island volumes, revealing that the ^-dependence 
of / reg is negligible compared to the ^-dependence of /hi er 
and /ch- We thus use / reg = p re g in the following and 
Eq. (Eh reduces to 

/ch + /hier = Pch ■ (6) 

We want to focus on the fraction /hi cr of hierarchical 
states, which we now determine from their influence on 
eigenfunction and level statistics. 
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FIG. 2. a) Successive magnifications of the phase space hierar- 
chy of the kicked rotor (K = 2.5). b) The fraction of hierarchical 
states /hier as a function of 1/h for K = 3.6,2.5,1.8 as deter- 
mined from eigenfunction (dots) and level statistics (diamonds). 
The solid lines show the predicted power laws according to Eq. (^|) 
and Fig. 2c. Data for K = 1.8 are extracted from table 2 in 
Rcf. c) Classical staying probability P(t) ~ t ' providing 

the values of 7 used in b) . The reason for the different fitting range 
for K = 1.8 is the longer lasting initial exponential decay. 

As a first quantity we determine 77, the density of an 
eigenstate averaged over a large area in the chaotic re- 
gion, far away from the islands. As long as this av- 
eraging area is large compared to Ti its specific choice 
and the used representation (e.g., Husimi) does not af- 
fect the following analysis. Normalization for 77 is chosen 
such that, an eigenstate uniformly covering the entire 
phase space would give rj = 1. The distribution of r\ 
(Fig. 1) shows a peak for small 77 stemming from the reg- 
ular states, intermediate values from hierarchical states 
and a Gaussian-like peak that we ascribe to the chaotic 
states. Assuming that the fraction / c h of chaotic states 
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extends on average uniformly (neglecting scar effects 1 19 
across a fraction / c h < p c h of phase space, we fit a Gaus- 
sian with area / c h and mean l// c h to the latter part 
of the 77-distribution. We find that the mean l// c h is 
clearly distinct from l/p c h (Fig. 1), the value one would 
obtain if there were eigenstates uniformly covering the 
entire chaotic region. The fit determines the fraction 
/hier — Pch — /ch of hierarchical eigenstates, which is 
shown as a function of h in Fig. 2b. 

Secondly, we want to quantify the influence of hierar- 
chical states on level statistics for mixed systems. To this 
end we extend the approach of Berry and Robnik poj , 
which assumes a random superposition of a regular (Pois- 
sonian) spectrum with measure pbr and a chaotic (GOE) 
spectrum with measure 1 — pbr : We take into account 
the hierarchical states and their spectrum. As most of 
them couple weakly to regular as well as chaotic states, 
their spectrum is superimposed independently to the rest 
of the spectrum. Typically, /hi er is smaller than / c h + / re g 
and therefore to first order these levels are randomly 
placed in the spectrum. Thus, effectively one can use the 
original Berry-Robnik approach, but one has to interpret 
the parameter /?br as the sum of regular and hierarchical 
fraction of states, /?br = /reg + /hier- We determine Pbr 
by fitting the nearest-neighbor level-spacing distribution, 
using the cumulative distribution as well as the so-called 
U-function fl2~l]| . The obtained values of /hi or from both 
fitting procedures agree within the error bars and their 
average as a function of h is shown in Fig. 2b. 

The two methods for determining /hior, eigenfunction 
and level statistics, rely on several assumptions, and one 
may not expect identical absolute values. Both meth- 
ods, however, show a clear power-law decay according 
to relation (||) for the fraction of hierarchical eigenstates 
as can be seen in Fig. 2b for three values of the kicking 
strength K, where 7 is extracted from the decay of the 
corresponding P(t) shown in Fig. 2c. 

We now want to derive Eq. (g). Let us remind that the 
classical power-law trapping P(t) ~ t f originates from 
partial transport barriers, e.g., Cantori, arranged in a hi- 
erarchical way in the chaotic part of phase space around 
islands of regular motion g. The simplest model ||] 
yielding P(t) ~ t~ 7 describes this hierarchy of the chaotic 
sea as a chain of volumes f2 n ~ u> n (u> < 1, n = 0, 1, . . .) 
connected by a flux $ n! „ + i ~ ip" (ip < 1) as shown in 
Fig. 3a. Then 7 = 1/(1 - In u/ In tp) holds §|. Quantum 
mechanically two neighboring volumes n and n + 1 are 
strongly coupled, as long as <E>„ in+ i > h, while they are 
weakly coupled in the opposite case 22 1 . This determines 
a critical flux $„*.„» + i = h with n* ~ ln?i/ Imp. The vol- 
umes £l n with n < n* correspond to the main part of the 
chaotic sea and support the chaotic eigenstates. Regions 
with n > n* correspond to the hierarchical part of the 
chaotic sea supporting the hierarchical states. Summa- 
tion of their volumes, finally, yields the fraction 
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FIG. 3. a) Sketch of the chain model with scaling volumes and 
fluxes, b) The fraction of hierarchical states /hier for three random 
matrix ensembles (7 = 1.15,1.3,1.86) as a function of l/h as de- 
termined from eigenfunction (dots) and level statistics (diamonds). 
The solid lines show the predicted power laws according to Eq. (0). 
c) Classical staying probability P(t) ~ t ' confirming the values 
of 7 used in b). 
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of hierarchical states. From 7 > 1 it follows, that in the 
limit h — > this fraction tends to zero, while the total 
number of hierarchical states goes to infinity. In case of a 
system with d 7^ 2 pairs of conjugate variables, one needs 
to replace Ti by h d ~ 1 in the above derivation. 

In order to check relation (]?]), we generated a random 
matrix model for the classical chain following the ap- 
proach of Ref. |22| . For each ensemble of these random 
matrices we fixed oj, ip and therefore also 7 and then var- 
ied h. We repeated the above eigenfunction and level 
statistics and extracted /hier (/reg — Prcg — in this 
model). The results for various 7 are shown in Fig. 3b. 
The numerical data agree very well with the expected 
power laws over three orders of magnitude in Ti and there- 
fore confirm relations (m) and (j^). 

In conclusion, we present numerical evidence that be- 
sides chaotic and regular eigenfunctions there exists a 
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third class, namely hierarchical eigenstates. Their frac- 
tion is determined from eigenfunction and level statistics 
and shows a power law as a function of h. We explain 
the origin of the power law and relate its exponent to the 



well known power-law trapping of mixed systems. 

We like to thank S. Tomsovic and L. Kaplan for helpful 
discussions. 
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